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We examine theoretically the effects of the bonding geometries at the gold-thiol interfaces on the 
inelastic tunneling spectra of propanedithiolate (PDF) molecules bridging gold electrodes and show 
that inelastic tunneling spectroscopy combined with theory can be used to determine these bonding 
geometries experimentally. With the help of density functional theory, we calculate the relaxed 
geometries and vibrational modes of extended molecules each consisting of one or two PDF molecules 
connecting two gold nanoclusters. We formulate a perturbative theory of inelastic tunneling through 
molecules bridging metal contacts in terms of elastic transmission amplitudes, and use this theory 
to calculate the inelastic tunneling spectra of the gold-PDF-gold extended molecules. We consider 
PDF molecules with both trans and gauche conformations bound to the gold clusters at top, bridge 
and hollow bonding sites. Comparing our results with the experimental data of Hihath et al. [Nano 
Lett. 8 , 1673 (2008)], we identify the most frequently realized conformation in the experiment as 
that of trans molecules top-site bonded to both electrodes. We find the switching from the 42 meV 
vibrational mode to the 46 meV mode observed in the experiment to be due to the transition of 
trans molecules from mixed top-bridge to pure top-site bonding geometries. Our results also indicate 
that gauche molecular conformations and hollow site bonding did not contribute significantly to the 
experimental inelastic tunneling spectra. For pairs of PDF molecules connecting the gold electrodes 
in parallel we find total elastic conductances close to twice those of single molecules bridging the 
contacts with similar bonding conformations and small splittings of the vibrational mode energies 
for the modes that are the most sensitive to the molecule-electrode bonding geometries. 

PACS numbers: 81.07.Nb, 72.10.Di, 73.63.Rt, 85.65.+h 



I. INTRODUCTION 

Electrical conduction through single organic molecules 
bridging a pair of metal electrodes has been studied ex- 
tensively, both experimentally and theoretically for more 
than a decade [1]. Knowledge of the microscopic de- 
tails of the bonding between the molecule and the elec- 
trodes is required in order to understand the transport 
properties of these systems. For molecules amine-linked 
to gold electrodes, it has been argued[2] that only one 
bonding geometry should be important. However, for 
the more commonly studied molecular wires with gold- 
thiolate bonding, it has been suggested that a variety of 
molecule-metal interface geometries may occur and influ- 
ence electrical conduction through these nanowires [1, 3— 
7]. Some possible geometries that have been proposed 
are those in which a sulfur atom of the molecule is lo- 
cated at a top site over a particular gold surface atom or 
over a bridge site between two atoms of the gold surface 
or over a hollow site of between three gold atoms. [8-11] 
However, which (if any) of these possibilities are actually 
realized in metal-molecule-metal nanowires has remained 
unclear [1]. Since a single molecule located between two 
electrodes and chemically bonded to both of them is not 
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accessible to scanning microprobes that are sensitive to 
atomic-scale structure, direct experimental determina- 
tion of the atomic scale geometries of the molecule-metal 
interfaces of metal-molecule-metal nanowires has been 
elusive and continues to be an important goal in this 
field [1]. In this paper, we show how inelastic tunneling 
spectroscopy (IETS) combined with ab initio computer 
simulations can distinguish between the different possible 
atomic-scale structures of the gold-thiolate interfaces of 
single-molecule nanowires and can therefore be used to 
definitively identify the molecule-electrode bonding ge- 
ometries that are realized experimentally. A brief ac- 
count of some of our preliminary results has already been 
published [12]. 

When the bias voltage applied across a molecular 
wire increases beyond the threshold value for the exci- 
tation of a vibrational quantum of the wire (to be re- 
ferred to as a phonon) at low temperatures, a conduc- 
tance step is observed in the current-voltage character- 
istic of the system. IEFS experiments have measured 
the molecular wire phonon energies in this way for var- 
ious molecules. [13-43] Theoretical studies of the IETS 
of molecular nanowires [12, 44-79] have also been car- 
ried out and have accounted for many aspects of the 
IETS data. These IETS experiments and comparison 
with the theoretical work have demonstrated conclusively 
that particular molecular species arc involved in electri- 
cal conduction through metal-molecule-metal junctions 
by exploiting the fact that the measured phonon ener- 
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gies provide a fingerprint of the molecule involved in the 
conduction process. The effects of changes in the molec- 
ular conformation [32, 41, 77] and in the orientation of 
molecules relative to the electrodes[59, 79] on the inelas- 
tic tunneling spectra have also been studied. However, 
the possibility that IETS might identify the bonding ge- 
ometries at the molecule-metal interfaces and thus re- 
solve the long standing "contact problem" of determining 
the atomic scale structure of the molecule-electrode inter- 
faces has only very recently begun to be explored. [12, 78]. 

We consider one of the simplest organic molecules, 1,3- 
propanedithiolate (PDT), bridging gold electrodes. Elec- 
trical conduction through this molecular nanowire and 
its IETS have recently been studied experimentally by 
Hihath et aZ. [29] using the statistical STM break junc- 
tion technique: By repeatedly forming a gold-PDT-gold 
molecular nanowire, measuring its conductance charac- 
teristic and then breaking the wire and reforming it, Hi- 
hath et al. [29] collected a large body of conductance and 
IETS data representing many separate realizations of this 
system. However, in common with previous break junc- 
tion experiments [1] , the specific geometry corresponding 
to any particular measurement was not identified. 

In the present study, we calculated the relaxed geome- 
tries, vibrational normal modes, zero bias conductances, 
and IETS of trans and gauche PDT molecules bonded to 
pairs of gold clusters in top, bridge, and hollow confor- 
mations. An important aspect of the methodology that 
we used is that we carried out relaxations of the entire 
extended molecules that we studied, relaxing the posi- 
tions of all of the atoms of the gold clusters as well as 
those of the atoms of the PDT molecule without any 
imposed constraints. This is necessary in order to ob- 
tain reliable results for the vibrational modes that have 
the largest amplitudes on the sulfur atoms that bond to 
the gold clusters and are therefore the most sensitive to 
the sulfur-gold bonding geometry. We also studied the 
dependence of our results on the sizes of the gold clus- 
ters and found good convergence of our results with in- 
creasing cluster size. Thus our results are expected to 
be relevant to PDT molecules bridging macroscopic gold 
electrodes as in experiments, and we do indeed find excel- 
lent agreement between our results and the experimental 
data of Hihath et al. [29] . We identify the most frequently 
realized conformation in the experiment [29] as that of 
trans molecules top-site bonded to both electrodes. We 
find the switching from the 42 meV vibrational mode 
to the 46 meV mode observed in the experiment [29] to 
be due to the transition of trans molecules from mixed 
top-bridge to pure top-site bonding geometries. We also 
conclude that gauche molecular conformations and hol- 
low site bonding did not contribute significantly to the 
experimental inelastic tunneling spectra [29]. For pairs 
of PDT molecules connecting the gold clusters in parallel 
we find total elastic conductances close to twice those of 
single molecules bridging the contacts with similar bond- 
ing conformations and small splittings of the vibrational 
mode energy of interest. 



The remainder of this paper is organized as follows: In 
Section II A, we outline how we calculate the equilibrium 
geometries and vibrational modes for the gold-PDT-gold 
molecular wires. In Section II B, we describe how the 
low bias elastic tunneling conductances of these molec- 
ular wires are calculated. In Section II C, we derive the 
perturbative expression that we use to calculate the con- 
ductance step heights associated with inelastic electron 
scattering in the molecular wires. In Section III A, we 
discuss the relaxed molecular wire geometries that we 
find. In Section MB, we present the results of our low 
bias elastic conductance calculations. In Section IIIC, 
we present the results of our calculations of the inelas- 
tic tunneling spectra of PDT molecules connecting gold 
electrodes and compare them with the experimental mea- 
surements of Hihath et al. [29] We present our conclusions 
in Section IV. 



II. THEORY 

A. Calculation of the low energy conformations 
and vibrational normal modes 

Our goal is to study inelastic tunneling processes that 
are sensitive to the structure of the gold-thiol inter- 
faces. These processes typically involve excitation of vi- 
brational modes with strong amplitudes on the sulfur 
atoms. Therefore, it is necessary to calculate accurate 
equilibrium geometries of the entire system being stud- 
ied, including both the molecule and the gold electrodes, 
and the phonon energies and atomic displacements from 
equilibrium for the vibrational normal modes of the whole 
system. We do this by performing our calculations for ex- 
tended molecules that consist of the molecule itself and 
two clusters of gold atoms (with up to 14 gold atoms per 
cluster) that represent the pair of electrodes to which the 
molecule binds. 

In order for this to be practical, we perform ab initio 
relaxations of these extended molecules based on density 
functional theory (DFT)[80, 81] with the GAUSSIAN 09 
[82] package. We obtain relaxed structures for the var- 
ious different bonding geometries and molecular confor- 
mations that we study by carrying out relaxations start- 
ing from different initial trial geometries. The relaxations 
are carried out with no constraints imposed on the geome- 
tries so that the relaxed geometries correspond to local 
minima of the total energy and therefore all of the com- 
puted vibrational modes about the relaxed geometries 
are well defined. 

We note that in most previous calculations of the in- 
elastic tunneling spectra of molecules bridging gold elec- 
trodes some or all of the gold atoms were frozen in the 
geometry of an ideal bulk gold crystal lattice. In that 
case the calculated geometries of the extended molecules 
were only partly relaxed and did not correspond to min- 
ima of the total energy. Therefore calculations of vibra- 
tional modes based on these partly relaxed geometries 
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suffered from a lack of rigor. Furthermore in the exper- 
iment of Hihath et al. [29] the molecular junctions were 
formed by crashing a gold STM tip into a gold substrate 
and then breaking the junction and repeating this pro- 
cess many times at cryogenic temperatures. This exper- 
imental "statistical STM break junction" methodology 
makes it very unlikely that the gold has the structure 
of an ideal bulk gold crystal anywhere near the PDT 
molecule that bridges the gold electrodes in the experi- 
ment. Thus models in which some of the gold atoms are 
frozen in the geometry of an ideal bulk gold crystal lat- 
tice are not expected to be appropriate for experiments 
of this type. This may explain why the authors of Ref. 78 
had difficulty modeling the results of IETS experiments 
on molecular junctions formed by the mechanically con- 
trolled break junction technique. 

In the present work, the atomic displacements from 
equilibrium in the normal modes and the corresponding 
frequencies and phonon energies are also calculated using 
DFT [82]. By carrying out systematic calculations for 
gold clusters of different sizes, we have checked that our 
conclusions are independent of the cluster size for the 
larger clusters that we study and thus are applicable to 
molecules bridging the nanoscale tips of experimentally 
realized macroscopic gold electrodes. 

We use the B3PW91 density functional (for the 
exchange-correlation energy) and the Lanl2DZ pseudo- 
potentials and basis sets. [82, 83] However, we have re- 
peated the calculations for some cases using the func- 
tional PBE0(PBE1PBE)[82, 84] for comparison and 
found little change in the result. We calculated the vibra- 
tional frequencies without using any scaling factors and 
will, in what follows, compare the corresponding calcu- 
lated phonon energies with the experimental data of Hi- 
hath et al. [29] without making any adjustments. 



B. Calculations of the elastic conductances in the 
limit of low bias 

Conduction through molecular nanowires in the low 
bias limit does not involve inelastic processes. We calcu- 
late the low bias conductances in this limit for the ab ini- 
tio extended molecule geometries obtained as described 
in Section II A using the extended Huckel model with the 
parameters of Ammeter et al. [85] to estimate the elec- 
tronic structures for these geometries. As is discussed 
in detail in Refs. 1 and 86, this methodology involves 
no fitting to any experimental data relating to trans- 
port in molecular wires or PDT molecules with gold con- 
tacts. Previous transport calculations based on extended 
Huckel theory have yielded elastic tunneling conduc- 
tances in agreement with experiment for other molecules 
thiol bonded to gold electrodes [86-90] and have also 
explained a number of transport phenomena observed 
in experiments on molecular arrays on silicon[91-93] 
and graphene nanoribbons[94] as well as electrolumi- 
nescence data, [95] current-voltage characteristics [95] and 



STM images[96] of molecules on complex substrates. Our 
calculated elastic tunneling conductances (at zero bias) 
for the gold-PDT-gold molecular nanowires are in good 
agreement with the values measured by Hihath et al. [29] 
We calculate the zero bias tunneling conductances g 
for gold-PDT-gold molecular wires from the Landauer 
formula 



9 = 9oJ2H({°}) 



(1) 



where go = 2e 2 /h and i|'({0}) is the elastic transmis- 
sion amplitude through the extended molecule at the 
Fermi energy Ep of bulk gold calculated within extended 
Huckel theory using the standard parameter set of Am- 
meter et al. [85] In the transmission amplitude, {0} means 
that the extended molecule is in a relaxed geometry cal- 
culated as is discussed in Sec. II A, i denotes an elec- 
tronic state of a carrier that is coming from the left lead, 
and j is the electronic state of a carrier that has been 
transmitted to the right lead. Vj and are the cor- 
responding electron velocities. The coupling of the ex- 
tended molecule to the electron reservoirs was treated 
as in previous work [86, 91-93, 97-101] by attaching a 
large number of semi-infinite quasi-one-dimensional ideal 
leads to the valence orbitals of the outer gold atoms of 
the extended molecule. This method of coupling the ex- 
tended molecule to the reservoirs is successful because 
most of the resistance of this system is due to the weakly 
transmitting molecule itself (as distinct from the cou- 
pling between the extended molecule and the ideal leads) 
and because the many ideal leads behave similarly to a 
large number of phase-randomizing Biittiker probes[102] 
in minimizing the influence of electronic dimensional res- 
onances due to the finite sizes of the gold clusters of the 
extended molecule. 

The transmission amplitudes tji were found by solving 
the Lippmann-Schwinger equation 



I*) = \$ ) + G (E)W\V) , 



(2) 



where |$o) is an electron eigenstate of an ideal semi- 
infinite one dimensional left lead that is decoupled from 
the extended molecule, Go(E) is the Green's function 
of the decoupled system of the ideal leads and the ex- 
tended molecule, W is the coupling between the extended 
molecule and leads, and |\&) is the scattering eigenstate 
of the complete coupled system. Since the basis set used 
in extended Huckel theory is non-orthogonal, the orthog- 
onalization procedure introduced in Ref. 103 was used to 
facilitate the above calculations. 



C. Perturbative theory of the IETS of molecular 



When the bias voltage V applied across a molecular 
wire exceeds tkj/e (where oj is the frequency of a vibra- 
tional mode of the system) , it becomes energetically pos- 
sible for electrons passing through the molecular wire to 
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emit phonons of that mode and a conductance step devel- 
ops in the low temperature current voltage characteristic 
of the wire. Since the bias voltages involved are small, 
the mode frequencies that we calculate at zero bias us- 
ing the ab initio approach described in Section II A are 
appropriate for estimating (via \V\ — huj/e) the values 
of the bias voltage at which conductance steps due to 
inelastic processes occur. 

In order to estimate the heights of these conduc- 
tance steps, we adopt a perturbative approach in the 
spirit of that proposed by Troisi et al.[A9) who trans- 
formed the problem of calculating IETS intensities into 
an elastic scattering problem. However, unlike Troisi et 
a?. [49], we formulate the present theory and results ex- 
plicitly in terms of elastic electron transmission ampli- 
tudes tj\ through the molecular wire. Such perturba- 
tive theories are appropriate at low bias voltages pro- 
vided that the elastic electron tunneling channel is well 
off resonance, [49] conditions that are satisfied for gold- 
PDT-gold molecular wires in the experiments of Hi- 
hath et al. [29] Since we are concerned with the case of 
weak (far off resonance) electron transmission through 
the molecule, inelastic forward scattering of electrons 
through the molecule is expected to make the dominant 
contribution to the steps in the conductance that are 
associated with the emission of phonons and the conduc- 
tance is expected to increase with increasing bias at these 
steps. [63, 70, 104] We shall assume this to be the case in 
the following derivation of the theory of IETS on which 
our numerical calculations of IETS intensities presented 
in Section III C are based. 

In the harmonic approximation the vibrational Hamil- 
tonian of an extended molecule is 



\p Pin 

^ 2m„ 



injmQinQjm > 



injm 



where Pi n and qi n are the i th atomic momentum and 
position coordinate operators for atom n. They obey the 
commutation rules [qimPjm] = ih~fii,jSn.m- The position 
coordinate q in is measured from the equilibrium position 
of the atom n. 

In terms of the creation and destruction operators 
aj, and a a for phonons of the normal modes a of the 
extended molecule with frequencies w„, the harmonic 
Hamiltonian takes the form 



H, 



vib 



and the atomic position coordinate operators are 



Qin — y ] \j - (di na a a + di na a) a ) 



2cj, 



(3) 



Here the coefficients d na are the displacements from 
their equilibrium positions of the atoms n of the ex- 
tended molecule in normal mode a normalized so that 
J2 n m nd* na , ■ d na — 5 a ' ia and m„ is the mass of atom n. 



Let us now approximate the interaction Hamiltonian 
U(r) describing the interaction between an electron at 
position r and the extended molecule by a sum of terms 



C/(r) = ^u„(r-R n ) 



(4) 



corresponding to the interactions of the electron with the 
individual atoms n of the extended molecule whose nuclei 
are located at R„. 

Then the change in the Hamiltonian U (r) if the atoms 
are displaced from their equilibrium positions R„ is 



AU(r) = u n (r - R„) - u n (r - R°). 



(5) 



To leading order in the displacements qi n of the atoms 
from their equilibrium positions, this becomes 



AU(r) = ^-q in — 



(6) 



Using Eq. (3 ), this yields the electron-phonon interac- 
tion Hamiltonian H eD 



H ep = AU(r) 



I h tV^j du„(r-R°) 

' 2^ V ^77 aa 2^ d ' na + hx - 



dri 



If A is a suitable scale factor and d na is the vector with 
components di na this becomes 

H ep = at,H a (r) + h.c. where 



H a (r) 



h 

2u> a A 



limV 

4->0 t— 1 



u n (r - (R° + Ad na )) - u n (r - R" ) 
A 



For the system of interest in this paper, the electron- 
phonon interaction H ep is a small perturbation. The op- 
erators a Q and a a induce transitions between states with 
differing numbers of phonon quanta while H a (r) simul- 
taneously gives rise to transmission of electrons through 
the molecular wire in the inelastic channel with scatter- 
ing amplitude . Assuming that the system is cryogenic 
so that the presence of any ambient phonons may be ig- 
nored, and that H a is so weak that t^ is effectively linear 
in the strength of H a , and that A is small enough that 



HJr) = 



h 

2^, 



\p u n (r - (R° + Ad na )) - «„(r - R°) 
^ A 



to a good approximation, it follows that 



ji ~ V 2w„ A ji ' 



(7) 



where t^ is the inelastic electron transmission amplitude 
that would be found if H a were 
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That is, is the change in the elastic electron 
transmission amplitude when the perturbation H a is 
added to the electronic Hamiltonian H equ n. static of the 
static molecular wire in its equilibrium geometry. But 
Hequii. static + H a is simply the electronic Hamiltonian of 
a static molecular wire whose nuclear positions R re are 
R n = R„ + Ad na , i.e., they are shifted from their equi- 
librium positions by vectors Ad na that are proportional 
to the displacements of the atoms from their equilibrium 
positions in normal mode a. 

This means that 



Relaxed Structures 



t a 



tf z ({Ad na })-tf t ({0}), 



where i|'({Ad nct }) is the elastic electron transmission 
amplitude through the molecular wire with each atom 
n displaced from its equilibrium position by Ad na where 
A is a small parameter and t|'({0}) is the elastic electron 
transmission amplitude through the molecular wire in its 
equilibrium geometry from state i with velocity Vi in the 
electron source to state j with velocity vj in the electron 
drain. It then follows from Eq. 7 that 



h 4({Ad na })-4({0}) 

lim — 

2uj a A^o A 




FIG. 1. (Color online.) Examples of relaxed structures of 
the trans-PDT molecule bonded to gold clusters in different 
ways. [105] (a): Both S atoms bond to the gold clusters in the 
top geometry, (b): One S atom bonds to a gold cluster in the 
top geometry while the other S bonds in the bridge geometry, 
(c): Both S atoms bond to the gold clusters in the bridge 
geometry, (d): Both S atoms bond to the gold clusters in the 
hollow geometry. Carbon, hydrogen, sulfur, and gold atoms 
are black, grey, yellow, and amber, respectively. 



calculations are carried out at the zero bias Fermi energy 
since the values of the bias at which the inelastic trans- 
mission occurs in the experiments of Hihath et al. [29] 
are low. 



The analog of Eq. (1) for the inelastic electron transmis- 
sion probability associated with the emission of a phonon 
of mode a is then 

T = y V A U I 2 = h lim y V A I ^({^W) ~ te Mi°}) I 2 

ij *i 

From this it follows immediately that the IETS intensity 
(i.e., the conductance step height) 8g a associated with 
the emission of phonons of mode a is 



5g a =go T a ,i.e., 
2e 2 h 

Sg a 



h 2uj a 
„2 



a-»o ^-r 1 Vi 



Vj |tjl({Ad„ a })-t|i({0})|2 



.or 



lim v Vj I te Mi Ad ^ - *M°>) | 2 

2-KWa A->0 Z—j v .\ A 



(8) 



Eq. (8) expresses the IETS intensities g a explicitly 
in terms of elastic electron transmission amplitudes t|j 
through the molecular wire. It states that, in the lead- 
ing order of perturbation theory, the scattering ampli- 
tude for inelastic transmission of an electron through the 
molecular wire is proportional to the change in the elastic 
amplitude for transmission through the wire if its atoms 
are displaced from their equilibrium positions as they are 
when vibrational mode a is excited. We evaluate if!- in 
Eq. (8) numerically to find the heights Sg a of the conduc- 
tance steps that arise from inelastic tunneling processes 
due to emission of phonons of vibrational mode a by ap- 
plying the methodology described in Section II B. These 



III. RESULTS 

A. Low energy conformations of the extended 
molecules 

As has been discussed in Section II A, we consid- 
ered extended molecules, i.e., molecules together with 
finite clusters of gold atoms to which the molecules 
bond, and relaxed these entire structures (the molecules 
and gold clusters) using density functional theory-based 
calculations. [80-84] No constraints were applied to the 
extended molecule geometries during the relaxation. 
This yielded calculated geometries that correspond to lo- 
cal energy minima of the extended molecules. 

Representative examples of the relaxed geometries that 
we obtained for extended molecules in which the PDT 
molecule has the trans conformation are shown in Fig. 1. 
Top, bridge and hollow-site geometries where each sul- 
fur atom of the melecule bonds to one, two, or three 
gold atoms, respectively, are shown in Fig. 1 (a), (c) and 
(d) . A mixed bridge-top structure where one sulfur atom 
bonds to two gold atoms and the other bonds to one gold 
atom is shown in Fig. 1 (b). 

For relaxed extended molecules having structures of 
the types depicted in Fig. 1, we found the extended 
molecules with the sulfur atoms binding to the gold clus- 
ters in bridge geometries to have the lowest energies. The 
hollow site bonding geometries had the highest energies. 
The top site bonding geometries had energies between 
those of the bridge and hollow cases. For example, the 
energy of the relaxed extended molecule in the bridge 
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geometry was 0.46 eV lower than that with the top ge- 
ometry and 1.33 eV lower than that with the hollow ge- 
ometry for our molecular wires with 13 Au atoms in each 
gold cluster. These numbers varied depending on the 
gold cluster size. But in all studied cases, the bridge 
geometry was more stable than the top geometry which 
was more stable than the hollow geometry for relaxed 
structures. 

We found that extended molecules for which DFT ge- 
ometry relaxations were started with the sulfur atoms 
over hollow sites on the surfaces of close packed gold clus- 
ters having atomic geometries resembling that of bulk fee 
gold invariably relaxed to bridge bonding site geometries. 
This is reasonable physically since the bridge site is not 
very far from the hollow site and a drastic rearrangement 
of the gold clusters is not required for a transition from 
hollow site bonding to a lower energy bridge site bonding 
geometry to occur. Thus while we were able to generate 
examples of relaxed extended molecule geometries (such 
as that in Fig. 1 (d)) with each sulfur atom bonding to 
three gold atoms (i.e., hollow site bonding) the structures 
of the gold clusters near these bonding sites were much 
more open than that of a defect-free surface of a fee gold 
crystal near a hollow site. Because of the much greater 
fragility and the higher energies of hollow site bonded 
structures relative to bridge and top site bonding, it is 
reasonable to expect hollow site bonding to be realized 
much less frequently than bridge and top site bonding 
(if it is realized at all) in statistical STM break junction 
experiments such as those of Hihath et al. [29] . As will be 
seen below, comparison of our theoretical inelastic tun- 
neling spectra with the experimental data of Hihath et 
al. [29] supports this expectation. 



B. Elastic conductances in the limit of low bias 

Our calculated zero bias conductance values for top, 
bridge, and hollow site bonded trans PDT molecules are 
in the ranges g =0.0012-0.0013g , 0.0015-0.0019g , and 
0.0032-0. 0095^0: respectively, for extended molecules at- 
tached to the larger gold clusters that we studied. Values 
for some specific gold cluster sizes are shown in Table I. 
All of these theoretical values agree with the experimen- 
tal value 0.006 ± 0.002 5o [29] to a degree that is typi- 
cal of the experimental and theoretical literature[106] for 
molecules thiol bonded to gold electrodes. We note that 
the calculations reported here did not involve any fitting 
to the data of Hihath et al. [29] or to any other molecular 
transport experiments. 

The calculated conductances show an increasing trend 
from the top to the bridge to the hollow bonding geome- 
tries due to the better electrical coupling between the 
molecule and electrodes for molecules bonding directly 
to larger numbers of electrode atoms. The calculated 
conductances for hollow site bonded molecules are closer 
to the experimental values [29] than are the calculated 
conductances for the the top and bridge site bonding. 



However, as we have already discussed in Section III A, 
structures with hollow site bonding are much more frag- 
ile and thus are less likely to be realized experimentally 
than bridge or top site bonded structures. Moreover, as 
will be shown below in Section III C, comparison between 
our theoretical predictions and the experimental inelas- 
tic tunneling spectra[29] provides strong evidence that 
top and bridge bonding were predominant in the experi- 
mentally realized systems of Hihath et al. [29] . Thus it is 
evident that the present elastic conductance calculations 
(in common with other elastic conductance calculations 
for molecules thiol bonded to gold electrodes that are 
available in the literature [106]) are not accurate enough 
for comparison between the theoretical and experimen- 
tal elastic conductance values to reveal which bonding 
geometries were realized experimentally. 

TABLE I. (Color online.) Calculated IETS intensities and 
phonon energies for Mode-I for one or two PDT molecules 
connecting gold clusters together with the calculated low bias 
conductances. T: top site bonding, B: bridge site bonding, 
and H: hollow site bonding of a molecule to a gold clusters. 
TT and BB indicate two molecules connecting the gold clus- 
ters in parallel and bonding in top and bridge geometries, re- 
spectively. Gauchel: One gauche C-C bond. Gauche2: Two 
gauche C-C bonds. Numbers in the configuration column in- 
dicate the numbers of atoms in the gold clusters. Notice that 
the low bias conductance decreases with increasing numbers 
of gauche bonds for similar bonding geometries, consistent 
with previous theories for longer molecular wires. [77, 107] 
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42.03 


9.54 


0.0016 


13H-H 13 Trans 


42.40 


34.02 


0.0095 


13H-T13Trans 


43.52 


11.74 


0.0029 


14B-B1 4 Gauchel 


39.06 


13.44 


0.00133 


14B-B14Gauche2 


36.58 


104.97 


0.00071 


14T-T14 Gauchel 


38.44 


1.80 


0.00034 


14T-T14 Gauche2 


36.58 


1.54 


0.00007 



a Experimental conductance measurement of Hihath et al. is 
about 0.006 g ±0.002 [29]. 



C. Inelastic Tunneling Spectroscopy of the 
Gold-Sulfur Interface 

In this Section we present the results of our calcula- 
tions of the inelastic tunneling spectra for trans PDT 
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molecules bonded to gold clusters in top, bridge, and 
hollow site conformations as well as for molecules with 
gauche bonds. We also consider some examples of pairs of 
molecules bridging the gold clusters. We compare our re- 
sults with the experimental data of Hihath et al. [29] and 
deduce which molecular and bonding conformations were 
most commonly realized in the experiment. In the cal- 
culations reported below, density functional theory was 
used to find the relaxed extended molecule geometries, 
vibrational modes and their frequencies as described in 
Sections II A and III A. The inelastic tunneling intensities 
Sg a corresponding to the various modes were calculated 
from Eq. (8), evaluating the elastic scattering amplitudes 
tjl that enter Eq. (8) as is described in Section II B. We 
will focus primarily on those vibrational modes that have 
the strongest amplitudes of vibration on sulfur atoms and 
the strongest IETS intensities since these modes are the 
most sensitive to the molecule electrode bonding geome- 
tries and are also prominent in the experimental IETS 
spectra. 



1. Top and bridge site bonding of trans PDT molecules 

Our calculations show that for trans PDT molecules 
the modes with strong amplitudes of vibration on the 
sulfur atoms that have the largest IETS intensities fall 
within the phonon energy range of a prominent feature 
of the experimental IETS phonon histogram [29] (shown 
in grey in Fig. 2(b)) that extends from 39 to 52 meV. 

The vibrational normal modes in this energy range and 
their phonon energies are shown in Fig. 2(a) for examples 
of extended PDT molecules with sulfur atoms bonded to 
14-atom gold clusters in top and bridge geometries. The 
vibrational modes are labelled I and II according to the 
nature of the atomic motion. The corresponding calcu- 
lated IETS spectra (IETS intensities vs. phonon energy) 
are shown in color in Fig. 2(b) for extended molecules 
with gold clusters of various sizes together with the ex- 
perimental IETS phonon mode histogram[29]. 

It is important to note that that vibrational modes I 
and II are strongly localized to the molecule itself and 
the gold atoms in its immediate vicinity. Their ampli- 
tudes of vibration decay rapidly in the gold clusters as 
the distance from the molecule increases. This is due 
to the fact that the frequencies of these modes are well 
above the phonon frequencies of bulk gold so that in the 
gold these modes are evanescent in character. For the 
largest gold clusters that we studied, the vibrational am- 
plitudes of modes I and II are very small already before 
the farthest gold atoms from the molecule are reached. 
Therefore, the calculated properties of these modes are 
not impacted significantly by the finite sizes of the gold 
clusters in the present study. 

The mode with the strongest calculated IETS inten- 
sities in Fig. 2(b) is Mode-I that is depicted in the top 
row of Fig. 2(a). In this mode (also known as the Au- 
S symmetric stretch mode) the sulfur atoms have the 
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FIG. 2. (Color online.) (a) Calculated vibrational modes 
in the phonon energy range from 38 to 53 meV for trans- 
PDT bridging gold nano-clusters with sulfur atoms bonded to 
gold in top-site and bridge-site geometries. Carbon, hydro- 
gen, sulfur, and gold atoms are black, grey, yellow, and am- 
ber, respectively. [105] Red arrows show un-normalized atomic 
displacements; the heavier arrows indicate the motion of the 
sulfur atoms. Mode-I has the stronger IETS intensities, (b) 
Calculated IETS intensities (colored) vs. calculated phonon 
energies for trans-PDT molecules linking pairs of gold clus- 
ters with between 11 and 14 Au atoms in each cluster. Re- 
sults are shown for both sulfur atoms bonding to the gold 
in top-site and bridge-site geometries and for top site bond- 
ing to one gold cluster and bridge site bonding to the other 
(bridge/top). The three ellipses enclose the calculated type I 
Mode- IETS spectra for pure bridge, pure top, and bridge/top 
bonding geometries for extended molecules with gold clusters 
containing various numbers of gold atoms. The experimental 
IETS phonon mode histogram of Hihath et al. [29] is shown 
in (darker, lighter) grey for (positive, negative) bias voltages. 
Modes of the types I and II in (a) are indicated by arrows. 



strongest vibrational amplitudes and move in antiphase, 
approximately along the axis of the molecule. Mode-II 
(also known as the Au-S antisymmetric stretch mode) 
is shown in the lower row of Fig. 2(a). It is similar to 
Mode-I except that in Mode-II the sulfur atoms move in 
phase with each other. As seen in Fig. 2(b) the calcu- 
lated IETS intensities for Mode-II are much weaker than 
those for Mode-I for both top and bridge site bonding. 
This difference between the IETS intensities of modes 



I and II can be understood physically as follows. Since 
in Mode-I the two sulfur atoms move in antiphase the 
gold-sulfur distances for both sulfur atoms either increase 
or decrease together as the extended molecule vibrates. 
These distances can be regarded as the widths of tun- 
nel barriers between the molecule and the two gold elec- 
trodes. The motions of the two sulfur atoms act in con- 
cert to widen or narrow both tunnel barriers together and 
thus to weaken or strengthen the electron transmission 
amplitude through the molecular wire. Therefore, the 
magnitude of the difference between the clastic electron 
transmission amplitudes of equilibrium and vibrating ge- 
ometries through the molecular wire in Eq. (8) is en- 
hanced. By contrast, in Mode-II, while the gold-sulfur 
distance for one sulfur atom increases, the gold-sulfur 
distance for the other sulfur atom decreases. These ef- 
fects of the motions of the two sulfur atoms on the elastic 
transmission amplitude through the molecular wire tend 
to cancel. As a result, the magnitude of the difference 
tj {({Adna}) — tjl({0}) that appears in Eq. (8) is smaller 
for Mode-II than for Mode-I and consequently the IETS 
intensity Sg a is seen in Fig. 2(b) to be much weaker for 
Mode-II. For this reason, although the amplitudes of the 
motion of the sulfur atoms in modes I and II are simi- 
lar and in both cases the motion is approximately along 
the molecular axis, Mode-I completely dominates the cal- 
culated IETS spectrum in the phonon energy range that 
corresponds to the main feature of the experimental IETS 
phonon histogram[29] in Fig. 2(b). 

That the IETS intensity for mode I is much larger 
than that for mode II can also be understood in terms 
of the Troisi-Ratner propensity rules [57, 58] for inelas- 
tic tunneling intensities. For a very simple model of a 
perfectly symmetric molecule with symmetric electrodes, 
Troisi and Ratner [58] found the IETS intensity for mode 

II to be exactly zero as a consequence of symmetry. The 
IETS intensity for mode II is not (in general) exactly zero 
in the present model because the electrodes are not sym- 
metric. We note that the electrodes are also not expected 
to be symmetric in current experimental molecular wire 
devices such as that of Hihath et al. [29] 

Notice that our theoretical results for phonon Mode-I 
of PDT molecules bonded to both electrodes in the top 
site geometry, bonded to both electrodes in the bridge ge- 
ometry, and bonded to one electrode in the bridge and to 
the other in the top geometry (the features around 45.5, 
40.5, and 42 meV in the theoretical spectra in Fig. 2(b), 
respectively) are reasonably well converged with respect 
to increasing gold cluster size. Also the Mode-I phonon 
energy ranges of the bridge site-bonded and bridge/top 
site-bonded wires do not overlap and are well separated 
from the phonon energy range of Mode-I for pure top site 
bonding. 

Based on these theoretical results, we propose the fol- 
lowing identification of the bonding geometries between 
the PDT molecules and gold electrodes that gave rise to 
the main features of the experimental IETS histogram[29] 
that is reproduced in Fig. 2(b): We identify the main 



peak in the experimental histogram [29] in Fig. 2(b) that 
is centred at ^46 meV as being due to trans-PDT 
molecules that bonded to both gold electrodes in the top 
site geometry. Thus we attribute the majority of the 
counts in the histogram to molecules in the pure top site 
bonding geometry. The weaker peak centred near 42 meV 
in the experimental histogram matches our predictions 
for molecules that bond to gold electrodes in the mixed 
bridge/top site geometry. The shoulder of the experi- 
mental histogram [29] at lower phonon energies centered 
near 40.5 meV corresponds to our results for molecules 
bonding to both electrodes in the bridge site geometry. 
Thus, our results show that IETS experiments can iden- 
tify experimental realizations of PDT molecular wires in 
bridge, bridge/top, and top site bonding geometries by 
measuring the energies of the phonons emitted during 
electron transport through these systems. 

Although the feature near 42 meV in the experimental 
histogram is not very prominent, Hihath et al. [29] ob- 
served a transition in which the IETS spectrum switched 
from exhibiting a peak near 42 meV to one near 46 meV 
as the molecular junction was stretched. They there- 
fore argued that the 42 meV and 46 meV vibrational 
modes correspond to distinct conformations of the molec- 
ular junction. However, they were not able to identify 
the conformations that give rise to these modes. Our 
theoretical results indicate that the switching from the 
42 meV to the 46 meV mode corresponds to the junction 
switching from trans PDT in a mixed top/bridge bonding 
geometry to trans PDT in a pure top bonding geometry 
as the junction is stretched. This is consistent with our 
finding that the pure top bonding geometry corresponds 
to a larger separation between the gold electrodes than 
does the top/bridge bonding geometry. 

In addition to modes I and II, the theoretical IETS 
spectrum in Fig. 2(b) shows other vibrational modes (la- 
belled "Mixed Modes") at lower energies, between 28 and 
36 meV. The calculated IETS intensities of these modes 
are weaker than those of Mode-I but in most cases not as 
weak as those of Mode-II. In the same energy range, the 
experimental histogram in Fig. 2(b) [29] shows very few 
counts. This suggests that the measurements of Hihath 
et al. [29] were not able to detect features in the inelastic 
tunneling spectrum having weak intensities and that the 
intensities of the "Mixed Modes" between 28 and 36 meV 
in Fig. 2(b) (as well as those of Mode-II) were at or below 
the detection threshold. 



2. Hollow site bonding 

For alkanedithiolates with longer chains of (6, 8, or 
10) carbon atoms bridging gold electrodes, the low 
bias conductances have been observed experimentally to 
drop by a factor of ~4 from the second lowest conduc- 
tance plateau to the lowest as the molecular junction 
is stretched. [5] It was conjectured [5] that this change 
in conductance may be due to the molecules switching 
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FIG. 3. (Color online.) Calculated IETS intensities (colored) 
vs. calculated phonon energies for irans-PDT molecules link- 
ing pairs of gold clusters with 12 (yellow) and 13 (green) Au 
atoms in each cluster for pure hollow site bonding and with 13 
Au atoms in each cluster for mixed hollow/top site bonding 
(magenta). The experimental IETS phonon mode histogram 
of Hihath et al. [29] is shown in (darker, lighter) grey for (pos- 
itive, negative) bias voltages. 



from hollow site bonding at one electrode and top site 
bonding at the other (hollow/top), to top site bonding 
at both gold electrodes. This conjecture is consistent 
with the expectation that the distance between gold elec- 
trodes should be larger and the conductance smaller for 
alkanedithiolate molecules bonded to gold electrodes in 
purely top site geometries than in hollow/top bonding 
geometries. However, whether it is correct or not has 
remained unclear in the absence of more direct experi- 
mental evidence. 

Can such a transition from hollow/top bonding to top 
site bonding plausibly account for the switch from the 
^42 meV phonon mode to the ~46 meV phonon mode 
that was observed experimentally [29] in the IETS spectra 
of the gold-PDT-gold system? As we have already men- 
tioned, in our DFT calculations PDT extended molecules 
started in hollow site bonded trial geometries in most 
cases relaxed immediately to bridge-site bonded geome- 
tries. This indicates that the hollow site bonding is much 
more fragile (and thus less likely to be realized experi- 
mentally) than bridge or top site bonding. 

Some examples of calculated IETS spectra for PDT 
molecules hollow site-bonded to one or both gold clusters 
are shown in Fig. 3. The calculated phonon energies and 
IETS intensities for these structures are more sensitive to 
the sizes of the gold clusters than they are for pure top, 
pure bridge or mixed bridge/top site bonding. As for 
the calculated elastic conductances of hollow site bonded 
structures (see Table I), this greater sensitivity is due to 
the more open structures of the gold clusters required for 
hollow site bonding and the resulting greater variability 
in the details of the hollow site bonding geometries. 

In Fig. 3 the energy range in which we find Mode-I 
phonons for pure hollow site bonding overlaps with the 



~42 meV feature of the experimental IETS histogram 
of Hihath et al. [29]. However, we find the structures 
with pure hollow site bonding to also support phonon 
modes with strong calculated IETS intensities (similar in 
strength to the Mode-I intensities for the top and bridge 
site bonded structures in Fig. 2(b)) in the mixed mode 
region below ~38 meV in Fig. 3 where the experimental 
IETS histogram shows few or no counts. This suggests 
that structures with pure hollow site bonding were re- 
alized rarely (if at all) in the experiment of Hihath et 
aZ.[29]. 

For mixed hollow/top bonding the calculated Mode-I 
phonon energy in Fig. 3 does not match either the ~42 
or the ~46 meV feature in the experimental histogram 
[29] but falls roughly half way between the two. Thus, in 
addition to the formation of bridge/top structures be- 
ing much more likely (based on our ab initio simula- 
tions) than the formation of hollow/top structures, the 
bridge/top bonding geometries account for the phonon 
energy of the 42 meV feature of the experimental [29] 
IETS spectra while the hollow/top bonding geometries 
do not. Thus it appears unlikely that hollow/top struc- 
tures played a significant role in the experiment of Hihath 
et a/. [29]. 



3. Gauche molecular conformations 

A gauche conformation of an alkanedithiolate molecule 
can be obtained starting from the trans conformation by 
rotating part of the molecule through ^60° about a C- 
C or C-S bond. Detection of gauche conformations of 
PDT molecules adsorbed on gold has been reported[108]. 
Also it has been conjectured [29] that the switching that 
Hihath et al. observed [29] in their inelastic tunneling 
spectra from the 42 meV vibrational mode the to 46 meV 
mode (as their gold-PDT-gold junction was stretched) 
might have been due to a change between gauche and 
trans configurations of the molecular junction. 

We have investigated theoretically how the presence of 
gauche bonds would affect the inelastic tunneling spectra 
of PDT molecules bridging gold contacts and the possi- 
ble implications for the experiment. [29] Representative 
results of our calculations are presented in Fig 4 (a) and 
(b) for PDT molecules with one and two gauche bonds 
connecting pairs of gold clusters in bridge site and top 
site bonding configurations. For the bridge site bond- 
ing configurations, the calculated energies of the vibra- 
tional modes of the gauche PDT molecules are not in the 
vicinity of either 42 or 46 meV. However, the calculated 
spectra exhibit vibrational modes with very strong IETS 
intensities at phonon energies between 32 and 38 meV, 
a range where the experimental histogram[29] in Fig. 4 
(b) shows very few counts. Thus it appears unlikely that 
these structures played a significant role in the experi- 
ment of Hihath et al. [29] . 

In Fig. 4 (b), the PDT molecule with one gauche bond 
and top site bonding does exhibit a vibrational mode 
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FIG. 4. (Color online.) (a) Representative examples of 
calculated sulfur-related vibrational Mode-I and II of gauche- 
PDT molecules connecting gold nanoclusters each with 14 
gold atoms, and the vibrational mode phonon energies. Gl 
and G2 label molecules with one and two gauche bonds respec- 
tively. The molecules connect to both gold clusters in pure 
bridge and pure top site bonding configurations. Carbon, hy- 
drogen, sulfur, and gold atoms are black, grey, yellow, and 
amber, respectively. Red arrows show un-normalized atomic 
displacements, (b) IETS intensities of gauche-PDT molecules 
linking pairs of gold clusters (14 Au atoms in each cluster) vs. 
the vibrational mode phonon energy. Sulfur atoms bond to 
the gold in top-top, bridge-bridge site geometries. The exper- 
imental IETS phonon mode histogram of Hihath et al. [29] 
is shown in (darker, lighter) grey for (positive, negative) bias 
voltages. Modes of the types I and II are indicated by arrows. 
* Scaled down by a factor of 10. 



with a phonon energy close to 46 meV. However, this is 



a type II mode with a low calculated IETS intensity. It 
is therefore unlikely to have been detected in the IETS 
experiment. [29] Other calculated vibrational modes of 
the same top-site bonded gauche structure in the phonon 
energy range from 27 to 36 meV in Fig. 4 (b)have some- 
what stronger but still weak IETS intensities and occur 
at energies at which the experimental histogram shows 
few if any counts. The PDT molecule with one gauche 
bond and top site bonding in Fig. 4 (b) also exhibits a 
vibrational mode with a phonon energy close to 46 meV 
but the calculated IETS intensity for this mode is much 
weaker still. These results indicate that gauche PDT 
molecules in top site bonding configurations also appear 
unlikely to have contributed significantly to the exper- 
imental IETS histogram of Hihath et al. [29]. We have 
also carried out calculations for gauche PDT molecules 
bonded to the gold electrodes in mixed top-bridge ge- 
ometries. The results of these calculations (not shown 
in Fig. 4) are qualitatively similar to those for the pure 
bridge bonding geometries of gauche molecules discussed 
above. Namely, these structures were also found not to 
have vibrational modes with calculated energies in the 
vicinity of either 42 or 46 meV but they were found to ex- 
hibit vibrational modes with very strong IETS intensities 
at phonon energies between 32 and 38 meV, where the 
experimental histogram[29] shows very few counts. Thus 
it appears unlikely that these mixed bridge-top bonded 
gauche molecules played a significant role in the experi- 
ment of Hihath et al. [29] 

The low bias elastic conductance measurements re- 
ported by Hihath et al. [29] also provide additional evi- 
dence indicating that switching between gauche and trans 
conformers is unlikely to have been responsible for the ex- 
perimentally observed transition [29] from the structure 
with the 42 meV vibrational mode to that with the to 
46 meV mode, as will be discussed in the following two 
paragraphs. 

We find the calculated sulfur-to-sulfur distances in 
PDT molecules with gauche bonds to be shorter than 
those for purely trans PDT molecules. Thus, if the tran- 
sition from the 42 meV vibrational mode to the 46 meV 
mode (that occurred as the junction was stretched) were 
due to switching between gauche and trans conforma- 
tions, then the 42 meV structure (being the shorter one) 
would correspond to the gauche geometry and the 46 
meV structure would have the trans geometry. However, 
the measured conductances of the PDT molecules ex- 
hibiting the 42 meV vibrational mode were larger than 
those of the PDT molecules exhibiting the 46 meV mode, 
as can be seen in Fig. 4 of Ref. 29. By contrast, re- 
cent theoretical work[77, 107] on other alkane-dithiolate 
molecules bridging gold electrodes found conductances of 
pure trans-molecules to be larger than those of molecules 
with gauche bonds, a finding consistent with the re- 
sults of our calculations of the low bias conductances for 
gauche and trans configurations of gold-PDT-gold molec- 
ular wires bonded to the electrodes in the same way; 
see Table I. These theoretical results suggest that if the 
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switch from the structures with the ~42 meV phonon to 
those with the ~46 meV phonon were due to a change 
from a gauche to a trans molecular configuration, then 
the switch would be accompanied by an increase of the 
low bias conductance of the junction, contrary to what 
was observed experimentally. [29] 

By contrast, our calculations show the low bias con- 
ductances of gold-trans-PDT-gold molecular wires in the 
bridge/top bonding configurations (with the ^42 meV 
phonon mode) to be higher than those for trans molecules 
in top bonding configurations (with the ^46 meV phonon 
mode); see Table I. Thus our interpretation of the ob- 
served transition from the ~42 meV phonon mode to 
the ~46 meV phonon mode being due to trans molecules 
switching from bridge/top bonding to pure top bonding 
is consistent with the experimental low bias conductance 
data[29]. 



4- Parallel Transport Through Pairs of Molecules 
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FIG. 5. (Color online.) The relaxed geometries and cal- 
culated conductances for pairs of trans PDT molecules con- 
necting gold clusters in parallel. 14 gold atoms per cluster. 
(a)-I Top site bonding. (a)-II Bridge site bonding. Calculated 
mode frequencies and IETS intensities are shown in (b) for 
the systems depicted in (a). The calculated mode frequencies 
and IETS intensities for single molecules similar to those in (a) 
connecting two gold clusters with the same kinds of molecule- 
gold bonding geometries are also shown for comparison. The 
experimental phonon mode histogram of Hihath et al. [29] for 
single molecules bridging the electrodes is also shown. The 
darker (lighter) grey represents the positive (negative) bias 
voltage. 



In the preceding discussion we have focussed on cases 
where a single molecule connects the two electrodes. 
However in statistical STM break junction experiments, 
such as that of Hihath et al. [29] , the electrodes can also 
be bridged simultaneously by two (or more) molecules 
in parallel. In this Section, we report theoretical results 
for the inelastic tunneling spectra of some examples of 
systems of this kind. 

In Figs. 5(a) I and II, we show the relaxed geome- 
tries for pairs of top site and bridge site bonded trans 
PDT molecules, respectively, connecting two gold clus- 
ters. As can be seen in Table I, in each case the calculated 
low bias elastic conductances of these pairs of molecules 
are roughly twice that of a single molecule connecting 
the contacts and bonded to them in a similar configura- 
tion. This means that the molecules in these systems are 
behaving as nearly independent parallel quantum con- 
ductors in a way reminiscent of parallel semiconductor 
quantum wires[109] connecting a pair of electron reser- 
voirs. Similar, roughly additive, behavior of the low bias 
conductances of other molecules bridging pairs of metal 
electrodes in parallel has also been found in previous the- 
oretical studies, for not too small intermolecular separa- 
tions [110-118]. 

The calculated inelastic tunneling spectra for these 
pairs of molecular wires are shown in Fig. 5(b). The in- 
elastic tunneling spectra for single molecules connecting 
the two electrodes with similar bonding geometries are 
also shown for comparison. The Mode I inelastic tun- 
neling spectra for the pairs of molecules are very similar 
to the corresponding spectra of the individual molecules 
apart from small phonon energy shifts and splittings, in- 
dicative of weak effective vibrational coupling between 
the molecules. However, for the pair of bridge site bonded 
molecules, the two Mode I features have noticeably dif- 
fering calculated IETS intensities. 

We also studied pairs of molecules with even smaller 
spatial separations than those of the molecules in Fig. 
5(a). Namely, we considered pairs of top site bonded 
trans PDT molecules binding to nearest neighbor atoms 
of both gold clusters, as distinct from the second neighbor 
gold atoms to which the molecules bind in Fig. 5(a)-I. 
As expected based on previous studies [110-118], the cal- 
culated low bias conductances of these systems differed 
more from twice the conductance of a single molecule con- 
necting the gold clusters than was the case for the more 
widely spaced molecules in Fig. 5(a). The splittings of 
the Mode-I phonon energies for these more closely spaced 
molecules were also considerably larger. However, we 
found that relaxed structures with such closely spaced 
molecules could only be generated for the smaller gold 
clusters that we studied; for the larger clusters DFT re- 
laxations resulted instead in structures similar to those 
in Fig. 5(a). This suggests that structures with such 
closely spaced pairs of molecules are unlikely to be real- 
ized in experiments such as those of Hihath et al. [29] . 

Experimental inelastic tunneling spectra for molecules 
bridging the contacts in parallel were not reported in 
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Ref. 29. Thus comparison between these theoretical re- 
sults and experiment, while of considerable interest, is 
not possible at present. 

IV. CONCLUSIONS 

We have calculated the relaxed geometries, zero bias 
conductances, and inelastic tunneling spectra of trans 
and gauche PDT molecules bonded to pairs of gold clus- 
ters in various ways and compared our results with the 
experimental STM break junction transport data of Hi- 
hath et al. [29] Although the gold clusters that represent 
the macroscopic electrodes in our calculations were not 
large (up to 14 gold atoms per cluster) they were suffi- 
ciently large for accurate modeling of the properties of 
the vibrational modes that we studied (as is discussed 
in Section IIIC 1) and our calculated inelastic tunneling 
spectra converged well with increasing size of the gold 
clusters. Thus, we expect our findings to be applicable 
to molecular wires bridging macroscopic gold electrodes. 
Consistent with this, our theoretical results agree well 
with the experimental data of Hihath et al. [29] We con- 
clude that IETS when combined with such calculations is 
able to reveal detailed, previously inaccessible, informa- 
tion about the atomic scale structures of the molecule- 
metal interfaces and thus resolve the long standing "con- 
tact problem" of single-molecule nanoelectronics. 

The success of our approach rests on the fact that 
ab initio density functional theory calculations of vibra- 
tional modes and their frequencies are believed to be rea- 
sonably accurate. The reason for this is that in the Born- 
Oppenheimer approximation they are electronic ground 
state total energy calculations and density functional the- 
ory has been designed and carefully optimized for calcu- 
lations of the total ground state energies of interacting 
electron systems. [1] While the results of transport cal- 
culations for molecular nanowires are well known to be 
less accurate, due in part to the well known difficulty 
of estimating the alignment of the molecular energy lev- 
els relative to the Fermi energy of the electrodes [1], we 
rely on transport calculations only for qualitative infor- 



mation such as the identification of the phonon mode in a 
particular frequency range that has the largest IETS in- 
tensity. Our identification of this mode is also supported 
by physical reasoning. 

We have definitively identified particular realizations of 
gold-propanedithiolate-gold molecular wires in a recent 
experiment [29] in which a trans PDT molecule bonded 
to a single gold atom of each electrode. We showed that 
this bonding geometry and molecular conformation was 
the one realized most frequently in the experiment. [29] 
Our theoretical results also showed that the switching 
from the 42 meV to the 46 meV mode observed in the 
experiment as the molecular wire was stretched [29] cor- 
responds to the junction switching from trans PDT in a 
mixed top/bridge bonding geometry to trans PDT in a 
pure top bonding geometry. 

Our results also showed that PDT molecular wires with 
gauche conformations connecting gold electrodes should 
have inelastic tunneling spectra that differ significantly 
from those with trans conformations. Comparing our re- 
sults with the inelastic transport data of Hihath et al. 
[29], we concluded that gauche conformations were real- 
ized at most rarely in their experiment [29]. 

We also found that if a pair of PDT molecules connects 
two gold electrodes in a parallel geometry the coupling 
between the molecules should affect the low bias elastic 
and inelastic transport characteristics of the system only 
weakly, unless the molecules are extremely close together, 
for example, if they bond in a top site geometry to adja- 
cent gold atoms of both electrodes. Structures with such 
closely spaced molecules were found to be unstable for 
the larger gold clusters in our DFT calculations. 
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